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I. INTRODUCTION 



Phenomenological modeling of interacting matter and (quantized) light in physics has a 
long and interesting history 1 . It is of quantum optical origin but is present in wide range 
of other branches of physics such as condensed matter ' or involving mechanical oscillators 4 . 
The Rabi model ', describing a qubit coupled to a single mode electromagnetic field, is the 
one which has attracted continuous attention for almost a century. More and less' recent 
studies on its integrability have inspired increasingly growing research. 

An existence of a symmetry of any quantum model is directly related to a quality of our 
understanding of its properties^. A 'sufficient' (in certain sense) symmetry can result in an 
integrability of the model . That is why seeking for any underlying symmetry of quantum 
models is always of great interest and often of great importance. In this paper we present 
our contribution to this activity. We consider a family of generalized single-mode Rabi 
models 10 : 



where a z and a x are the Pauli matrices, a and u correspond to the energy gap of the spin 
and boson, respectively, whereas a and are the annihilation and creation operators of 
quantized mode of light satisfying canonical commutation relation, [a, a'] = I. It is assumed 
that the coupling between the qubit and the field, controlled by the strength constant g, 
incorporates k > photons. 

In this paper, by solving an operator Riccati equation associated with Eq. (1), we con- 
struct an operator exhibiting significant similarities to the parity operator acting on the 
bosonic space. This operator, the generalized parity, can be used to simplify multi-photon 
Rabi model (1) and transform it to a block-diagonal form. Our work is a complemen- 
tary expansion of certain results obtained in Ref. 10 for k = 1 and k = 2 in the context of 
approximate methods of solving the Rabi model. 

The paper is organized as follows: In Sec. II we present operator Riccati equation asso- 
ciated with (1) serving as a main tool applied in our studies. Next, in Sec. Ill the known 
results concerning the k — 1, 2 cases are reviewed. Sec. IV has been devoted to the con- 
struction of the generalized parity and contains main results of our work. Finally, in Sec. V, 
followed by conclusions, we apply the general construction to a simple example. 




(1) 
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II. A TOOL: RICCATI EQUATION 



Multi-photon Rabi model considered here belongs to a general class of qubit-environment 
composite systems described by Hamiltonian 



H QE = H Q ® I E + I Q <g> H E + H 



hit 



H + V 
V 1 " H 



H 



QE) 



(2) 



where Hq (He) is the Hamiltonian of the system (environment). H int is the interaction of 
the qubit with its surroundings. Iq and Ie are identities acting on corresponding Hilbert 
spaces C 2 and He- The total Hamiltonian Hq E acts on C 2 <g> He and the symbol ~ should 
be understand as Ht corresponds to 1 in the sense of block operator matrix representation of 
operators. This correspondence is established via the isomorphism C 2 <8> "He ~ "He © He- 
Finally, the form of remaining operators H± and V depends upon how Hq, He and H int are 
defined. 

Any steps toward diagonalization of Hqe is valuable as it can be followed by variety of 
different approximation schemes ". There is often an additional benefit emerging form such 
transformations which can help to exhibit useful symmetry properties being often obscured 
by an 'improper choice' of a basis. As it is pointed out below it is also the case of the 
multi-photon Rabi model (1) discussed in this paper. 

Our idea originates from an observation that Hamiltonian Hqe can be converted to a 
block-diagonal form 



s 1 h qe s 



H + + VX 
H 





- (VX)t 



with S 



Ie -X f 



provided that X satisfies an operator Riccati equation 



XVX + XH + - H_X - V f = 0. 



(3) 



(4) 



For general considerations regarding an operator Riccati equation we refer the reader to 11 ' 12 . 
This equation provides valuable tool allowing to study the exact diagonalization , stationary 
states 11 and in general, the dynamics 15 of two level open quantum systems 16 . From the 
decomposition (3) it is evident that the dynamics of a qubit-environment quantum system 
is actually governed by the Riccati (4) and pair of uncoupled Schrodinger equations. 
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For the ^-photon Rabi model studied in our paper 

H± = ua ] a ± (g*a k + g(a^) k ) , V = aI Hs (5) 

and the corresponding Riccati equation reads as follows: 

aX 2 + XH + -H_X-a = 0. (6) 

Its mathematical properties has already been addressed in literature 1 '. In Eq. (6), H± are 
operators acting on the bosonic Fock space "Hb, a is a real constant, whereas X is a solution 
to be found. If it does not lead to a confusion, we write a rather than aIn B , with being 
the identity on Hb- 

III. KNOWN SOLUTIONS: k = 1,2. 

For the sake of self-consistency, we begin with reviewing known solutions and their prop- 
erties for the two particular cases, where k = 1, 2. For the simplest possible case, k = 1 the 
solution of the Riccati equation (6) was found in 13 to be the bosonic parity operator 

P = 5>*>>(n| = 5>l)»H (7) 

ngN n€N 

which can also be written in a more compact form as P = exp(^7ra^a), where {|n)}„, g N is 
the Fock basis, i.e., a)a\n) = n\n). Such operator is both hermitian and unitary, hence it is 
an involution (P 2 = I^ B ). Interestingly, it solves Eq. (6) for both a = (dephasing 18 ) and 
a^fl (exchange energy between the systems is present) cases, although they reflect quite 
different physical processes. 

In the context of RWA-type approximation methods the two-photon (k = 2) Rabi model 
was studied in details within 10 . The two-photon parity operator 

(8) 

was introduced therein. It has not been stated explicitly in but the parity T is, as will be 
shown below, a solution of the Riccati equation (6) for k = 2. 

IV. GENERAL CASE: k > 



T := exp 



i— a)a (a) a — l) 



. 2 



In what follows we show how to construct a solution of the Riccati equation (6) with 
coefficients H ± provided by (5) in the general case k > 0. Before we start let us emphasize 
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that the parity operator P (T) introduced in the preceding section solves Eq. (6) not only 
for k = 1 (k = 2) but also for all odd k = 2n + 1 (even, of the form k = In + 4) cases. This 
has already been noticed in 10 . Here we will not only fill the remaining gap k = 2n + 2 but 
also present unified approach allowing to obtain a linear solution for arbitrary k. As a first 
step toward constructing this solution, we define a family of orthogonal projectors 

oo 

P; := \n,l)(n,l\, with |n, /) := \kn + 1 - 1), (9) 

n=0 

for n G N and 1 < I < k. The states \n, I) satisfy the following orthogonality condition: 

(^|j;0 = $kn+i—l,km+j-l = fiij^nrn, (10) 

where 8 xy is the Kronecker delta. The first equality in Eq. (10) comes from the orthogonality 
of the Fock basis. The second one can be justified as follows. 

When i = j both sides of (10) reduce to S nm since 5k n +i-i,km+i-i = 8 nm . If i ^ j (say 
i > j) the right hand side is zero. The left hand side also vanishes, as one gets either m = n 
or m ^ n in this case. Indeed, if to = n then, to get nonzero left hand side, one would 
expect i — j = what is impossible. Finally, for m ^ n (m > n, say) one would expect 
that k(m — n) = j — i (in order to keep the left side nonzero) what also does not occur as 
k > i — j and it cannot divide % — j. 

For a given family of orthogonal projectors one can split the space Hb into k subspaces 
so that 

k 

n B = n 1 ®n 2 ®---®n k - 1 ®n k = ^ni, (11) 

i=i 

where "H/ := P; ("Hr)- The symbol © indicates the (orthogonal) direct sum of Hilbert spaces. 
Hereafter, we use it interchangeably with + when it refers to the sum of operators. 

The decomposition (11) allows us to think of operators H± as of k x k block operator 
matrices [H^j^xfe such that 

R l ^:=P l (R ± )P m :'H l ^n m . (12) 

Obviously, H" act within the space Hi, and they may be considered as operators H± re- 
stricted to the space Hi, i.e., H" := (H±)m r The off-diagonal elements H^ act between 
the subspaces Hi, H m and therefore transform state from one space into the other. In this 
simple picture instead of speaking of compositions, sums and any other operations involving 
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two or more operators acting on T-Lb, we operate with corresponding k x k matrices. The 
block operator matrix representation of H± is useful provided that it has a relatively simple 
form. As an example, an occurrence of a block diagonal structure (H^ = i( m Hf) would be 
an ideal situation. As it will be shortly seen this is indeed the case here. 

Since (10) holds, it is a matter of straightforward calculations to show that (i, l\a^a\j, m) = 
^jmSij6i m and (i, l\a k \j, m) = r n jm 5 i j- 1 5i m , for a given k, where 



(kj + m — l)! 

Z jm = kj + m-l and rj jm = ^ ^ - ^ + - _ ^ . (13) 
By making use of the above equations we obtain 

oo / oo \ 

P;(a f a)P m = Y {(h l W a \j, m )) \i,l)(j,™>\ = ^2Cjm\jJ)(j,m\ 5 im 

ij=o \j=o J (14) 

= (km + (I - l)I Hl ) 8 lm = N t Si m , 
where we have introduced N/ - the number operator restricted to the subspace Hi, : = 



{a)a)\ Hl : 



N, = fcn, + (J-l)I«„ m = ^n\n,l)(nX (15) 

n=0 

with Ifi t being the identity on Hi. In a similar fashion we have 

oo / oo \ 

Pi(a k )P m = ((hl\a k \j,m)) \i,l)(j,m\ = Ij^Vjmlj ~ l,l){j,m\ I S lm , (16) 

i,i=0 \j=0 ) 

from which Ai := J2'jLo r lji\j ~ 1 , (j, ^| is nothing but (a k )\u l and plays on Tii a role of 
annihilation operator. 

Combining all the above results into a single equation we finally obtain the block matrix 
representation of H ± with respect to the decomposition (11). It has a diagonal structure 
indeed: 



wNz ± (g*A l + gAj)] 5 lm = E^6 lm , (17) 



for l,m < k. 

Having (17) in place, we can formulate the result of this paper: 

Proposition 1. The solution to the Riccati equation (6) is given by the generalized parity 
operator 

k oo 

X fc = ^^(-l)>,Z)<n,Z|, for k>0. (18) 



1=1 n=0 
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Proof. We begin by denning partial parities J; to be 

oo 

j r - e i™ l = J2(-l) n \ n ,l)(n,l\, l<k. (19) 

n=0 

It is justified to refer to J; as a parity since it possesses all the desired properties required 
from the parity operator on Hi, namely 

J? = [N,, h] = 0, and J Z A;J ; = -A,. (20) 

As an immediate consequence of this conditions JjH® J/ = H^, what in the block operator 
matrix terminology developed in the preceding section leads to 

X fc H + X fc ~ diagfJiH^Ji, . . . , JfcHfjfc] ~ H_, (21) 

or simply X fc H + = H_X fc . Therefore, to prove (6) holds true, it is sufficient to show that X k 
is an involution. It can be established in a following way: 

k 

3^~diag[J? l ... J Ja~0I« I =I« B . (22) 

i=i 

Both in Eqs. (21) and (22) we have taken into account the correspondence X fc ~ diag[Ji, . . . , □ 
V. EXAMPLES 

It is interesting to see how the solutions which have been found in the previous section 
can be used to recover the known results for k = 1,2. When k = 1 there is only one subspace 
of *Hb, namely Hb itself, and Xi is equal to the bosonic parity operator (7). For k = 2 there 
are only two projection within the family of operators (9), i.e., 

oo oo 

Pi = ^|2n)(2n|, P 2 = J2 \2n + l)(2n + 1|, (23) 

n=0 n=0 

or in a compact form Pi^ = |(I-h b ^ P)> which according to (11) split the bosonic Hilbert 
space into two subspaces. The first one consists only of odd, while the second one of even 
Fock states: 

Kb = span{|2n) : n G N} © span{|2n + 1) : n E N}. (24) 
The block operator matrix representation of H± with respect to (24) reads 

~ H T , (25) 



H± ~ 



Hi 1 ) 







Ji 














Ji 





Hi 2) 







h 







Hf 




J 2 
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where diagonal entries are explicitly given by (17). To see that the two-photon parity T is in- 
deed the same as X2 one only needs to invoke a simple fact, (— l)"^ 2 ™ -1 ) = (_ i)™( 2n + 1 ) = (— 1)™, 
then 

00 00 
T = ^(-1)^^)71) (n| = ^(-l) n (|2n)(2n| + |2n + l)(2n + 1|) . (26) 

n=0 n=0 

Eq. (25) also allows us to identify X 2 ~ diag[Ji, J 2 ]. 

VI. CONCLUSION 

In summary, we have found the solution of the operator Riccati equation (6) associated 
with the multi-photon Rabi model (1) indicating certain symmetry of the original model. 
This solution is a natural candidate for a parity operator as it not only reduces to well 
known one- 13 and two-photon JI) parities but also its properties are that of a typical parity 
operator. 

We have also proved by an explicit construction that the solution of (6) with coefficients 
given by (5) exist for every k > 0. We have excluded from our analysis the case where 
k = only because it does not reflect any relevant physical phenomena. Since there is no 
interaction between the systems they evolve in time separately. 

The symmetry associated with the generalized parity Eq. (18) allows to convert Rabi 
model into a block-diagonal form. Such a transformation decouples the original qubit-boson 
eigenproblem into a pair of bosonic Schrodinger equations. Despite certain mathematical 
subtleties of the multi-photon Rabi models 19 we hope that the results presented in this 
paper can serve as a starting point for useful approximation methods. 

Finally, let us notice that any involution J for which JH + = H_J is a solution of (6) for 
k > 0. Unfortunately, the question whether all the solution to (6) are of that special kind 
still remains open. 
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